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Abstract 

The quantum Fourier transform has been implemented on a three bit nuclear magnetic resonance 
(NMR) quantum computer, providing a first step towards the realization of Shor's factoring and 
other quantum algorithms. Implementation of the QFT is presented with fidelity measures, and 
state tomography. Experimentally realizing the QFT is a clear demonstration of NMR's ability to 
control quantum systems. 

PACS numbers 03.67.-a, 03.67.Lx, 02.70.-c, 89.70.+C 



Quantum computers are devices that process 
information in a way that preserves quantum co- 
herence. Unlike a classical bit, a quantum bit, 
or 'qubit,' can be in a superposition of and 
1 at once. This nonclassical feature of quan- 
tum information allows quantum computers to 
perform some computations faster than classical 
computers. For example, quantum computers, 
if constructed, could factor large numbers more 
rapidly [||, search data basis more quickly |||], 
and simulate quantum systems more efficiently 
^] than is possible using current classical algo- 

rithms i i i i i i Q. 

A key subroutine of algorithms for factor- 



ing and simulation 0] |T^] IQ is the quan- 
tum Fourier transform (QFT) |]l5| jl^ |17]. In 
essence the QFT takes a 'position' state \x) to 
the corresponding 'momentum' state \p) and is 
defined as follows: 
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Where q is the dimension of the systems Hilbert 
space. 

In general the QFTq transforms the input am- 
plitudes as, 

QFT,Y.f{x)\x)^Y.hv)\p). (2) 
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Where the coefficients f{p) are 



To implement the QFT, these gates, 
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For example, the two qubit QFT corresponds to 
the unitary operator. 
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This operator shows the QFT separating the in- 
put states by degrees in the first row and col- 
umn, and then by 90 degrees, 180 degrees and 
270 degrees, multiples of ^. 

Equation (4) shows that the QFT has effects 
similar to that of the classical Fourier transform. 
In particular, if /(a) is periodic with period r, 
then /(c) will exhibit a spike at c = r. This is the 
key to Shor's algorithm which allows a quantum 
computer to factor very large numbers in poly- 
nomial time. The classical Fourier transform re- 
veals the periodicity in functions, the QFT re- 
veals periodicity of wavefunctions. 

As formulated by Coppersmith, the QFT can 
be constructed from two basic unitary opera- 
tions, Aj, operating on the jth qubit 
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and Bjk operating on the jth and kth qubits 
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are implemented on the lead bit, j = L — 1. Re- 
peating the above sequence of gates to all L bits 
as j is indexed from L — 1 to will complete the 
QFT. This sequence of quantum logic gates can 
be realized NMR. The idea of using nuclear spins 
as the basic unit of a quantum computer was 



proposed by Lloyd [|18[ , and detailed schemes for 
using NMR as a method of quantum computing 



were proposed by Cory et al ||l^ and Gershen- 
feld and Chuang |^0[. In NMR a series of ra- 
dio frequency pulses are used to control the ex- 
cess magnetization of an ensemble of quantum 
states. NMR experiments are easily visualized 
by picturing the excess magnetization as a vec- 
tor pointing in some direction and the pulses as 
rotations about the various axes. In addition, a 
bilinear coupling term in the Hamiltonian allows 
for quantum superposition. 

The Hamiltonian of a three spin (qubit) NMR 
sample with J-coupling is 



H 
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where /j = crj/2. The three bit QFT was im- 
plemented via NMR using the three carbon-13 
spins of an alanine sample. The resonant fre- 
quency of carbon-13 at 9.4 Tesla is approxi- 
mately 100.617MIIz. The carbonyl was labeled 
spin 1, Ca was labeled spin 2, and C/3 spin 3. The 
chemical shift of the three alanine carbons are 
12587IIz, OHz, and -3435IIz respectively. Cou- 
pling constants between the three spins are J12 
= 54Hz, J23 = 35Hz, and J13 = 1.2Hz. Relax- 
ation time Ti for alanine is approximately 1.56s 
while T2 is about 420ms. 

The Aj matrix described above can be broken 
up into idempotents E^ — E^ + ax[Ej^ + E^). 



The pulse sequence of the Aj gate can now be de- 
termined using the geometric algebra formalism 

Hi, 
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This pulse program reads: apply a pulse along 
the y-axis that rotates spin j 90 degrees, apply a 
pulse along the x-axis that rotates j 180 degrees. 
Magnetization on the z-axis would be rotated to 
the positive x-axis. Since this experiment starts 
with the spins at thermal equilibrium (pointing 
along the z-axis) the above sequence for the Aj 
gate can be replaced by the simpler ^ pulse along 
the positive y-axis. 

The Bjk gate, which can be constructed using 
the coupling between qubits, In terms of idempo- 
tents reduces to l-E-^E-'^+e^'^E-^E-'^. Again 
using geometric algebra this yields the following 
pulse sequence: 



i-y, - (2^) - i-)i 
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The notation 0/27r Jj^ represents an interval of 
spin evolution under the coupling Hamiltonian. 
The final three pulses effectively perform a ro- 
tation around the z-axis. These pulses are not 
necessary, however, since the same effect may be 
achieved by rotating the prior pulses of the ex- 
periment. 

The complete pulse program is the combina- 
tion of Aj and Bjk gates described above. In 
this implementation, the necessity of performing 
a swap gate has been removed by reordering the 
bits at the appropriate interval. 

The complete pulse program is. 
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(lefe) - W-x - (si) - (^)x 
(at) - (^)^ - (^)-x - (^)-x - 

W-.-(^)'-(f)!-(^)-.- 



This sequence includes a number of (vr) pulses 
to refocus couplings during the intervals they 
should be inactive. 

The pulse sequence takes advantage of knowl- 
edge of the starting state of the system at the 
beginning and end of the program by replacing 
Hadamard transforms with ^ pulses. In the mid- 
dle of the sequence the full Hadamard was indeed 
used. 

Figure 1 shows selected theoretical and exper- 
imental spectra following the quantum Fourier 
transform of the state l\ + II + /f on the three 
qubit NMR quantum computer. 

The fidelity of the QFT calculated using the 
measure 
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is 87%. Here p is the density matrix minus 
the part that is proportional to the identity (in 
NMR, this is called the 'reduced' density ma- 
trix; it should not be confused with the reduced 
density matrix got by partially tracing the den- 
sity matrix for a composite quantum system over 
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Figure 1: The three carbon-13 atoms of alanine 
after performance of the QFT. The top spectra 
are theoretical while the bottom are experimen- 
tal. Peaks in NMR spectra show the difference 
in energy level of single spin flips. Each spin has 
four peaks since its energy level is dependent on 
whether the other two spins are up (along the 
magnet) or down. This shows up clearly in spin 
B which has resolved J-couples to both of the 
other spins. The J-coupling between the A and 
C spins is very small and, therefore, the four 
peaks are not totally resolved. These peaks tell 
the magnitude of only some of the terms of the 
density matrix. 



some of its subsystems). This measure reflects 
both imperfections in the applied pulses and de- 
lays, as well as decoherence. To a first approx- 
imation, decoherence during the course of the 
QFT attenuates the entire density matrix. This 
is shown in figure 2. Therefore, we can approxi- 
mately separate the errors caused by experimen- 
tal imperfections by renormalizing p^xp to its at- 
tenuated average. Using this the fidelity of the 
operations themselves is above 98% over the 6 
gates in (11). 

The fidelity of 87% corresponds to an error 
rate of 97.7% over the six gates which, while 
high, does not attain the error rate of 10~^ 



required for robust quantum computation |22]. 
These errors arise primarily from spatial inho- 
mogeneities in the radio frequency fields which 
we believe can be improved. 

In conclusion, using NMR, the QFT has been 
implemented on a three bit quantum system and 
the fidelity with which we can transform an ini- 
tially diagonal state has been measured. Al- 
though the fidelity does not reach that required 
for fault tolerant computing, it is easily high 
enough to permit studies on small quantum sys- 
tems including quantum simulations. A particu- 
larly straightforward use of the QFT is in quan- 
tum chaos: as Balazs and Voros |Q pointed out, 
a simple version of the quantum baker's map 
can be performed by QFTs and Schack [^] has 
shown how such a quantum map might be real- 



ized on a quantum computer [25|. 

The authors thank S. S. Somaroo and C. H. 
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